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1 Introduction 



Quark matter at high density provides an interesting laboratory for the controlled study of 
non-perturbative physics [1J. Novel phenomena such as color superconductivity and color- 
flavor locking have been shown to result from physics near the quark Fermi surface. Rigorous 
results can be obtained in the limit of asymptotic density. Less, however, is known about 
what happens at intermediate densities, where the effective QCD coupling is still large. The 
intermediate density region is more likely to be realized in neutron star cores or laboratory 
experiments. 

The 't Hooft anomaly matching conditions || constrain the realizations of chiral sym- 
metry in the low energy phase of a gauge theory. In this paper we investigate whether these 
conditions, orginally derived at zero density, can be used to constrain the behavior of matter 
at non-zero density. There are several obstacles to this generalization. A basic observation 
is that Lorentz symmetry, which is broken at non-zero matter density, plays an important 
role in the original derivation of these matching conditions 0, [3|, |J, and in the physics of 
the anomaly. In gauge theories at zero density, unbroken global chiral symmetries imply the 
existence of massless spin-1/2 states. As elaborated in [§, £|, ||, these degrees of freedom are 
responsible for the IR singularities associated with the anomaly. It is not at all obvious that 
the same considerations apply to the quasiparticles comprising the low-energy spectrum at 
finite density. Yet, as noted by Sannino ||, the various phases of color superconductivity do 
indeed satisfy the anomaly matching conditions. 

't Hooft's original argument for his anomaly matching conditions involves the use of 
"spectator" fermions to cancel any anomalies which result from the gauging of flavor sym- 
metries. In the low energy limit of the theory, one is left only with the massless states 
(composite or elementary) and the spectator fermions. If this long-wavelength theory is to 
be consistent, the anomaly generated by the spectators must be cancelled by that of the 
other massless states. This argument can be repeated in the case of finite density, with gap- 
less quasiparticles playing the role of massless fermions in the low-energy effective theory. 
However, due to the existence of a Fermi surface there is a large degeneracy of zero energy 
states. Naively, one might expect an enhancement of the anomaly prefactor by an amount 
proportional to the area of the Fermi surface, which scales like /i 2 (the square of the chemical 
potential). Clearly, we need to understand more precisely how quasiparticles contribute to 
the anomaly. 

We will show that the anomaly matching conditions (AMCs) continue to apply at finite 
density, and constrain the possible quasiparticle spectrum |J. In essence, the singularity 
structure of the three-current correlator implied by the anomaly must be reproduced by the 
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effects of quasiparticles in the finite density theory, thereby constraining their flavor quantum 
numbers. Applying the Landau-Cutkowsky rules, we see that only quasiparticles satisfying 
special kinematic conditions contribute to the singularity - there is no degeneracy factor due 
to the Fermi surface. 

We focus our attention on cold but dense matter theories while a discussion of the anomaly 
at finite temperature can be found in Ref. 0. 

This letter is organized as follows. In section 2 we discuss the computation of the anomaly 
at finite density and show that it is unaffected by the presence of a chemical potential. 
Further, the anomaly continues to imply the existence of IR singularities, even at non-zero 
density. In section 3 we use the Landau-Cutkowsky rules to show that these singularities (in 
the absence of spontaneous symmetry breaking) require the existence of gapless quasiparticles 
in the low energy spectrum. 



2 Anomaly at Finite Density 



In this section we show that the anomaly is unaffected by the presence of a chemical potential, 
as are the implications for singularities of the three current correlator. Heuristically, density 
effects are IR in nature whereas the anomaly can be computed from the UV behavior of the 
theory. Points (A) and (B) below follow simply from this observation, although (C) does 
not. 

A. The standard derivation of the anomaly from the point of view of UV divergencies 
involves the careful treatment of the variation of the fermionic partition function || 

dM(i>,i>) exp i J d A xijj(i J^jijj , (1) 

under the infinitesimal change of variables 

tp — > (1 + ia(x)'j 5 ) ip , 

7p — ► i>(l + ia(x)^) . (2) 

Although the action in the exponent of the integrand is invariant under (g), the fermionic 
measure dM is not: dM. — > J~ 2 dM. The Jacobian determinant J is evaluated using 
gauge invariant regularization in the basis of the eigenstates of Dirac operator p, i.e. ip(x) = 
J2n a n(fin(x) , 4>(x) = En«n</>nO) and dM = Uda m dd m . 

Note, that if the chemical potential is not zero, one may still define the measure in 
the basis of the eigenstates of Dirac operator with zero chemical potential. All further 
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considerations follow the same path and produce the familiar result 

J = exp[-ij d 4 xa(x) ( ^e^F^F pa (x) ) ] (3) 

This implies that the anomaly equation for the axial current will take the same form at 
non-zero fermion density: 

<W = - (4) 

B. It is known that the anomalies arise via quantum corrections, (i.e. renormalization). 
In particular axial anomalies are associated, in a diagrammatic language, with triangular 
diagrams. 

The three point function in electrodynamics is: (the generalization to non-Abelian theo- 
ries is straightforward) 

T^x {h, h, q) = i J d A Xl d A x 2 (0|T (V^(x 1 )V u (x 2 )A x (0)) | 0) e ik ^ +ik ^ , (5) 

with q = ki + k 2 and V and A vector and axial currents respectively. Since the anomaly 
is independent of a fermion mass, we set it to be zero. The classical axial Ward identity 
<l X T pv x = is modified according to (for a more detailed discussion see M before and after 
formulae (6.31-32)): 

q% uX = Ag> + Ag> , (6) 



with 



A« = / ^Tr 



(2tt)< 



(7) 



and A^) is obtained by interchanging [i — ► v as well as k\ — > k 2 . Note that if we could shift 
the integration variable p to p + k 2 in the second term of ((?]) the AW term would vanish 
identically. Since the integrals are linearly divergent, it can be shown that (at zero density) a 
translation of the integration variable produces extra finite terms ruining the classical Ward 
identity. 

It is natural to ask what happens to the triangle anomaly at finite matter density. In 
principle, the contributions from triangle diagrams could depend on the chemical potential 

We compute the difference A^(/j) — Aj^(0) parameterizing a possible deviation in the 
anomaly calculations. We explicitly evaluate the finite density effects for the i — 1 case (the 
effects for the i = 2 are identical): 



(2tt)< 



^57^-7 — 7, 7m - -7757^-7 — 7r7/x- 



i> - ^70 i> - h - ^70 i> i> - h 



By applying a momentum shift (p to p — u), with the shift vector u T = (/i, 0, 0, 0), we expect 
a possible non null contribution to be encoded in the following integral (see ||): 



8A$M = u r J^-jjpfbMM* 0) 

and 

ff u u \ -a \ V P (g ~ k iY (P~ hY (p~ qY ~\ nn , 

f[p, fa, k 2 ) = -i Ae paflv — —2 - — — -2 . (10) 

lp z (p-ki) [p — K2)(p — q) 
Applying Gauss's theorem for the case of four- dimensional Minkowski space, we have: 

5A$(ja) = u T J — i —f(p,k 1 ,k 2 )=u T j^lim o p 2 p T f(p,k u k 2 ) 

= ^W^l^p^-^ + j (11) 

= ^{Wkl + u°k%) = 0, 

p T pP g T P 

and, in the last step, we used lim — — = — . 

p^oo 4 

Thus, for any finite density we have no extra contribution to the axial anomaly with the 
respect to the zero density one. 

Adler and Bardeen [|K| showed, at zero density, that the axial anomaly coefficient is one 
loop exact. We expect this theorem to hold at finite density, although we will not attempt 
to provide a proof. Heurisitcally we can say that the introduction of a density term does not 
affect the superficial degree of divergence of the higher-order triangle diagrams with respect 
to the zero density case. Since the latter have a lower degree of divergence with the respect to 
the simple triangle diagram we have no momentum-routing ambiguity and hence vanishing 
contribution to the anomaly. 

C. At zero density the anomaly implies an IR singularity in the three-point function ([5]) 
involving vector and axial currents 0. 

Now the form of T^ v \ is strongly restricted by Lorentz, Bose and permutation symmetries: 



T apfl (fa, k 2 , q) = + F x (k{ - k T 2 ) e Tam + F 2 k la k"k^e a/ 3 P ^ - k 2p k°k 2 e a0 , 



(J/./ 



+ ^3 



(Jfi 



+ 



(12) 



In general ( |T2| ) may contain tensor-like structures (denoted here by . . .) which unlike pseudo- 
tensorial ones do not contribute to the anomaly. 
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At non-zero matter density T pu \ contains more terms, due to the presence of a constant 
Lorentz four- vector u = (/i, 0). These terms will, however, be still restricted by the rotational 
50(3), Bose and permutation symmetries. It is convenient to introduce the 4- vector rj with 
u = fir]. Due to the presence of a new independent 4-vector rj the F functions, are general 
functions, i.e. F = F(ki-rj, k 2 -r), fi, k\, k 2 ). Following Ref. ]7) the most general decomposition 
in terms of invariant amplitudes is: 



(h,k 2 ,q) 



+Fi (k\ — k 2 ) € Tapp + F 2 ki a kik 2 e a/ 3 Pfl 



k 2a k-y k 2 ( 



afipfi 



k\ p k-y k 2 ( 



a[3ap, 



+ F 4 



V, 



tv 2p tv 1 hv 2 t a ^ ap 

rk-^ k 2 ( 



afipp 



T) p k 1 k 2 ( 



a[3ap 



+ F 5 k la k^r] l3 e a /3 P p + k 2p klrfe afi(Jll + F e k 2a kyT] p e a p pll + k Xp k 2 rfe 



+ F 7 



rjak^Sappn + V P k 2 V P ^af3ap 



+ F 8 



afBap 

kipk^f] 6 a f3 ap + k 2a k 2 rf € a f3 pp 



+ F 9 [k 2p k^r]^e af3ap + k la k^e a/3pp: j + F w [r^pk^e^^ + -q^rf ] t aPpp 
+ [Fa {kipkia - k 2p k 2a ) + F 12 {k Xp q a - r) p k 2a ) 



+F 13 {k 2p 7] a - r] p k la )} k^k 2 rfe a 



13-yp 



(13) 



One can explicitly verify that the previous amplitude is even under the simultaneous ex- 
change of p — > a and k\ — > k 2 . In the previous expression we dropped the flavor structure of 
the three point function since it is not relevant for our discussion. 

In what follows we fix the kinematics such that k\ = k\ = k 2 with k\ • r\ = k 2 ■ r\ = ko- 

Further constraints arise when imposing the anomaly equation along with the conserva- 
tion of the vector currents (i.e. in formulae): 



q' x T apfl (ki,k 2 ,q) 
klT^^k^q) 



2 2 1 fi/ 2 fc O/3crp 5 

k 2 T app (kt, k 2 ,q) = . 



(14) 
(15) 



After imposing these constraints and conveniently dotting the final results in the vectors 
rj, k\ and k 2 (see Ref. 0) we have the following relations: 



2Fj + (F 5 + F 6 - F s - F 9 ) k + (F r - F w ) 

(F 5 - F 9 ) k 2 + (F 6 - F 8 ) h ■ k 2 + (F 7 - F 10 ) k 
(F 5 -F 9 )k 1 -k 2 + (F 6 -F 8 )k 2 + (F 7 -F 10 )k = 0. 



2tt 2 







(16) 



The latter equations are due to the anomaly condition in Eq. (|l~5"l). The vector current 
conservation, on the other side, leads to: 



Fx + F 3 h ■ k 2 + F 2 k 2 + F 4 k 







F 6 h ■ k 2 + F 5 k 2 + F 7 k = , 
F 8 h ■ k 2 + F 9 k 2 + F w k = . 



;i7) 
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Combining the previous relations we have the following anomaly constraint: 



2F 3 k 1 ■ k 2 + 2k 2 F 2 + 2k F 4 - -±r*-F 



2 k 2tt 2 ' 

with F = F 5 — F 9 = F 6 — F 8 . 

In the absence of Lorentz breaking (i.e. zero density) F 4 = F = and a singularty in 
q 2 = (k\ + k 2 ) 2 — > arises when considering the limit kf = k 2 = k 2 = 0: 

lim F 3 = -4-? . (19) 
fc 2 ^o 3 2n 2 q 2 K ' 

This pole is interpreted as due to the pion when chiral symmetry is spontaneously broken. 

What happens at finite density ? Following Ref. [0] let us define the following null vectors 
ki = fco(l,0, 0, 1) and k 2 = k (l,sm8,0,cos8). In these variables we have: 

F 3 = ! *± FO + cos^) 

3 47T 2 fc2(l-COS0) fc o (l-COS0) (1-COS0) ' 1 ' 

Let COS6 1 7^ 1, and consider the double limit k — > and q 2 — * 0. We deduce (g 2 = 
2A; 2 (1 -cos0)): 

lim F 3 = . (21) 

fc 2_ 0;go _ 2 7T 2 g 2 V ; 

The form of this singularity is the same as the zero density singularity, as long as we consider 
the ko — > limit. We again interpret this singularity as due to massless excitations. In the 
case of unbroken chiral symmetry it implies the existence of gapless excitations with the 
charges necessary to match the zero density anomaly. 

It is important to note that this singularity does not imply that the anomalous decay 
7r — > 27 is unaffected by finite density effects (see 0). In medium, the gauge boson acquires a 
plasma mass, and hence the physical decay takes place at non-zero energy. At such kinematic 
points we cannot isolate the singularity cleanly, and density effects are possible. 

In order to better clarify how the infrared singularities emerge at finite density, we in- 
vestigate the singularity properties of quasiparticle scattering in the following section. 



3 Quasiparticles and Singularities 

Since the anomaly equation is unmodified by finite density effects, we can apply 't 
Hooft's spectator fermion construction to conclude that quasiparticles must generate the 
same anomalies as the fundamental fermions. In this section we examine in more detail how 
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this occurs. We are particularly interested in why the Fermi surface degeneracy does not 
alter the result, and why massive excitations do not contribute to the anomaly. 

As discussed above, the anomaly at finite density still implies a 1/Q 2 pole in the function 
Fi{Q 2 ). We can use the Landau- Cutkowsky rules [p| , to identify the intermediate states 
which are capable of providing this singularity. We first prove that the singularity structure 
is not modified by matter density effects. Then we investigate the kinematics associated 
with the singularity to gain a better understanding of why this is the case. 



Using the Landau-Cutkowsky rules [|TTJ , the discontinuity in the triangle graph is given 
dpo d 3 p 



by: 

Disc T^ x oc 



(2.) {P + h ~ 

e (p° + k° 1 -^)5[{p + k 



p (p + h + k 2 - u) a (p - u) T (p° - fijS (p- u) 2 

e (p° + fc° + k\ -fi)s [<j> + h + k 2 - uf] , 

(22) 



u 



times the trace factor Tr ^^y^y^y . Here we used the result [II] that the leading 
singularity in the physical region is given by the graph(s) in which each internal particle is on- 
shell and propagates forward in time. In other words, we replaced each internal propagator 
by 27ri6(Zo — /J>)S (h — p) 2 — i , where I is the particle momentum. 

We can now perform the integral in po by using the first 5 function, which requires that 
Po — fi = \p\. We obtain 

d 3 p 



Disc T^x oc J (P + k i) P (P + h + k 2 ) a (p) T 

(\p\ + fc?) 5 [(p + hf] (\p\ + k° 1 +k 



(p + ki + k 2 ) 



(23) 



with the 4- momentum p = (\p\,p). This integral no longer has any explicit dependence on 
the finite density term fi, and is identical to the integral obtained in the zero density case. 
We see that p — fi in the finite density case merely plays the role of the energy p Q in the 
zero density case. Hence the singularity structure is unchanged at finite density. 

Finally, we would like to understand why the Fermi surface degeneracy does not affect the 
anomaly calculation. To this end, let us examine the on-shell kinematic conditions satisfied 
by the quasiparticles in (^). Suppose we regard the triangle graph as a spacetime process, 
with the k\ vertex earliest in time. Let the incoming gauge boson momentum k\ be null. In 
order that the quasiparticles with momenta p and p + k\ both be on-shell, p and k\ must 
be parallel. Thus, a particular point on the Fermi surface is selected by k\. the rest of the 
surface cannot contribute to the singularity (see figure [[]). (One of the particles emerging 
from the k\ vertex is actually a quasihole on the opposite side of the Fermi surface from 
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the quasiparticle. Which is which depends on whether k\ and p are aligned or anti- aligned.) 
Finally, we note that if we take the energy of the gauge boson k\ — > 0, we will not be able 
to produce quasiparticles whose energy spectrum has a gap: only gapless quasiparticles can 
reproduce the anomaly in the entire physical region. 

A 

P • 




P II k 

1 

Figure 1: The particle- hole pair contributing to the anomaly singularity is determined by 
the external momentum k±. 

4 Discussion 

As shown, 't Hooft's anomaly matching conditions are as valid at non-zero matter density 
as at zero density. They provide a non-trivial check on recent results on QCD at asymptotic 
density. Morever, they apply even at strong coupling, where dynamical calculations cannot 
be reliably performed. 
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